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Abstract. We wish to assess the relationship between the population of planetary 
nebulae (PNe) and a given parent stellar population from a theoretical point of view. 
Our results rely on original population synthesis models used to estimate the expected 
luminosity-specific PN density accounting for different evolutionary scenarios and star 
formation histories, as observed in galaxies in the near Universe. For a complete PN 
sample, we find that 1 PN/1.510 6 Lq is a safe (IMF-independent) lower limit to the 
traced global bolometric luminosity of the parent stellar population. A tentative appli- 
cation to Virgo cluster data allows us to place a lower limit at ~ 7% for the global B 
luminosity of the cluster provided by "loose" intergalactic stars. 



1 Introduction 

Recent works on diffuse light in nearby clusters have renewed interest in plan- 
etary nebulae (PNe) as tracers of the parent stellar population. The efficient 
detection of these nebulae in the field, even at 70 kpc or more beyond the optical 
radius of bright galaxies (see e.g. Peng et. al 2004), makes them ideal candidates 
to trace the "loose" stellar population possibly filling the "empty" space among 
galaxies. PNe can be detected up to in nearby clusters at a distance of 20 Mpc, 
like for instance for the Virgo or the Fornax cluster (see Arnaboldi et al. 2003; 
Feldmeier et al. 2003). 

A major issue in this regard concerns the link between PN number density 
and global (bolometric) luminosity of their parent stellar population. This ra- 
tio (the so-called "luminosity-specific PN density" or, shortly, the "a ratio" , as 
defined first by Jacoby 1980) provides in fact direct information on the distinc- 
tive physical properties of the underlying stellar population as far as its PN 
component can be picked up by the observations. 

Population synthesis models can be especially useful in this sense, allowing 
us to tackle the problem from a fully theoretical point of view, thus overcoming 
most of the natural limits of the empirical calibrations (i.e. observational bias, 
incomplete sampling etc.). On this line, we want to give here a brief progress 
report of our project aimed at calibrating "a" for different evolutionary scenarios 
and star formation histories as observed in the local galaxy sample and likely 
consistent with the expected properties of the diffuse stellar component of the 
intergalactic medium. 
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Fig. 1. Theoretical fuel consumption for stars along the Hot-PAGB evolution according 
to the models of Paczynski (1971; solid curve) and Vassiliadis and Wood (1994; solid 
dots and triangles). Fuel is expressed in Hydrogen-equivalent solar mass. The stellar 
core- mass range to produce PNe is marked, after Dorman et al. (1993). 



2 Theoretical fundamentals 

Theory of simple stellar populations (SSPs), as developed in Renzini & Buzzoni 
(1986) and Buzzoni (1989) (hereafter RB86 and B89, respectively), provides the 
basic reference framework for our task, allowing us to evaluate the value of "a" 
for coeval and chemically homogenous stellar aggregates. 

Briefly, the expected number (A^pn) of PNe for a SSP of total (bolometric) 
luminosity L to t is: 

Ap N = BL tot r PN . (1) 

In the equation, B is the so-called "specific evolutionary flux" (a nearly age- and 
IMF-independent parameter of the order of 210~ n [L^yr" 1 ], see B89), and 
tpn is the PN lifetime (i.e. the time for the nebula to be detectable in [O in] 
and/or Ha surveys). 

The PN lifetime directly modulates with the stellar core-mass evolution, the 
latter depending on the nuclear fuel burnt at high temperature (T er r — 10 5 K) 
along the Hot Post-AGB (H-PAGB) phase (this is a direct result of the so- 
called "Fuel Consumption Theorem" of RB86). Figure 1 summarizes the model 
predictions, after Paczynski (1971) and Vassiliadis and Wood (1994), displaying 
the H-PAGB fuel consumption for stars of different core masses. 1 In the plot, 
we express the fuel in terms of Hydrogen-equivalent solar mass (i.e. 1 g of H- 
equivalent mass = 6 10 18 ergs; cf. RB86). 

The inferred H-PAGB lifetime for SSPs of different age, according to B89, is 
reported in Fig. 2. It can be defined as thpagb — Fuel/^n-PAGB, being -^h-pagb 

1 Note that actual PAGB core mass of stars can sensibly be lower than the original 
stellar mass when leaving the Main Sequence Turn Off point. This is a consequence 
of the mass loss processes that heavily affect stellar evolution along the red and 
asymptotic giant branch phases. 
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Fig. 2. Hot-PAGB core mass lifetime, for SSP evolution according to B89 synthesis 
models. A Salpeter IMF is assumed, and a Reimers mass loss rate coefficient r\ — 0.3. 
The explored metallicity spans the range from Z ~ 1/20 — > 2 Zq, as labeled in the 
plot. The PN dynamical timescale (Ta yn ) is also marked, for reference. 
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Fig. 3. Theoretical luminosity-specific PN density for SSPs of different metallicity. 
Models are according to B89, with the distinctive parameters as in Fig. 2. 



the luminosity of the stellar core at the onset of the nebula ejection (see e.g. 
Paczynski 1971 for details). It is important to note that, in general, thpagb is 
shorter than the dynamical timescale for nebula evaporation (Td yn — 30 000 yr; 
cf. Peimbert 1990); this implies that, at least for SSPs of young and intermediate 
age, tpn should actually coincide with thpagb- Only for issp <; 12 Gyr we have 
that tpn — Tdyn, so that the luminosity-specific PN density could eventually be 
written as 

a = — ^- = Bt pn = B inf{ thpagb, T dy „} (2) 
-kssp 

From eq. (2), and taking the B89 synthesis code as a theoretical reference, 
in Fig. 3 we report the expected luminosity-specific PN density for SSPs of 
different metallicity. Some physical age constraints exist for SSPs to produce 
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Fig. 4. Expected luminosity-specific PN density for template galaxy models (labeled 
with the corresponding Hubble type), according to Buzzoni (2002). Each curve tracks 
galaxy evolution from 5 Gyr to 15 Gyr (small solid dots), as labeled for the case of 
Ellipticals. Big triangles report the empirical estimate of a = Npn/L sa i from the PN 
census of Jacoby (1980) for the galaxies in the Local Group. The case of the Milky 
Way is singled out in the plot by the big solid dot. 



PNe; from one hand, for tgsp <J 10 8 yr SSP evolution is dominated by stars with 
M* ;> 5 M Q ending up as Supernovae. On the other hand, no PNe are expected 
to form from stars with actual mass M» < 0.52 M (Dorman et al. 1993; see 
Fig. 1); this is likely to occur at some point along SSP evolution for ages between 
10 <J tssp <, 18 Gyr, the exact details depending on the efficiency of the mass 
loss processes 2 . 



3 The luminosity-specific PN density in the galaxies: 
implications for the intergalactic medium 

The theoretical estimate of the luminosity-specific PN density is important when 
we infer total bolometric luminosity of their parent stellar population in the 
intergalactic medium. So far, we have relied upon the empirical calibration of a 
on the bulge of M31, with all the uncertainties induced by the poorly understood 
effects of metallicity, age and star formation history. 

Our new results, based on this improved theoretical approach, indicate that 
a luminosity-specific PN density 

a = 1 PN/1.510 6 L (3) 

2 The B89 SSP models parameterize mass loss via the coefficient 77 according to the 
Reimers (1975) standard notation. A value of r\ = 0.3 is adopted here, as an indicative 
value for the observed mass loss in the Galactic globular clusters (see B89 for details) . 
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is a safe estimate for old SSPs. Quite interestingly, this value docs not strongly 
depend on metallicity and, even more importantly, on the stellar IMF (as £>, in 
eq. (2), is nearly insensitive to the IMF details, as shown by RB86). 

On the basis of the derived SSP evolution, and matching the Buzzoni (2002) 
template galaxy models, we have further elaborated our synthesis scheme com- 
puting in Fig. 4 the expected evolution of a for the different star formation 
histories of galaxies along the Hubble morphological sequence. A comparison 
with the observed PN population for the Local Group members (Jacoby 1980, 
including the relevant case of the Milky Way) is also reported in the figure. 

As a tentative application of our results, we could match the recent obser- 
vations of the Virgo cluster by Arnaboldi et al. (2003) in order to assess the 
luminosity contribution from the intergalactic stellar background. Some 37 PN 
candidates have been detected by Arnaboldi et al. (~ 25% of which possibly 
consisting of spurious Lya emitters at z = 3.14, as claimed by the authors) 
across an area of 0.196 deg 2 on the cluster. As the bona fide PNe only span the 
brightest magnitude bin of the standard luminosity function we had to correct 
for sample incompleteness according to Ciardullo et al. (1989). 

By adopting the reference range of a for 15 Gyr galaxy models of Fig. 4, an 
intergalactic surface brightness (in bolometric) could be derived as 26.0 <; ^boi <, 
28.5 mag arcsec -2 , depending on the assumed star formation history for the PN 
parent stellar population. 3 This implies a B surface brightness for intergalactic 
stars in the range 27.3 <; <, 30.3, mag arcsec -2 . Compared to the smooth 
B surface brightness distribution of the bright-galaxy component of Virgo (i.e. 
^b — 27.4 mag arcsec -2 , see Arnaboldi et al. 2003), this allows us to place a 
confident lower limit to 7% for the cluster global luminosity to be provided by 
"loose" intergalactic stars. 

References 

1. M. Arnaboldi et al.: Astron. Journal 125, 514 (2003) 

2. A. Buzzoni: Astrophys. Journal Suppl. Ser. 71, 817 (1989) [B89] 

3. A. Buzzoni: Astron. Journal 123, 1188 (2002) 

4. R. Ciardullo, G. H. Jacoby, H. C. Ford, J. D. Neill: Astrophys. Journal 339, 53 

5. B. Dorman, R. T. Rood, R. W. O'Connell: Astrophys. Journal 419, 596 (1993) 

6. J. J. Feldmeier, R. Ciardullo, G. H. Jacoby, P. R. Durrell: Astrophys. Journal Suppl. 
Ser. 145, 65 (2003) 

7. G. Jacoby: ApJ 42, 1 (1980) 

8. B. Paczynski: Acta Astr. 21, 47 (1971) 

9. M. Peimbert: Rev. Mexicana Astron. Astrof. 20, 119 (1990) 

10. E. W. Peng, H. C. Ford, K. C. Freeman: Astrophys. Journal 602, 685 (2004) 

11. D. Reimers: Mem. Soc. Roy. Sci. Liege, 6th Ser. 8, 87 (1975) 

12. A. Renzini A. Buzzoni: 'Global properties of stellar populations and spectral evo- 
lution of galaxies'. In Spectral Evolution of Galaxies, ed. C. Chiosi, A. Renzini 
(Reidel, Dordrecht 1986) pp. 195-235 [RB86] 

13. E. Vassiliadis, P. R. Wood: Astrophys. Journal Suppl. Ser. 92, 125 (1994) 



3 For our calculations we adopt a distance modulus (m — M) = 30.8 mag for the Virgo 
cluster, and an absolute bolometric magnitude M bo i,0 = 4.72 for the Sun. 



